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Abstract 

The origin of the power-law decay measured in the power spectra of low Prandtl number Rayleigh-Benard 
convection near the onset of chaos is addressed using long time numerical simulations of the three-dimensional 
Boussinesq equations in cylindrical domains. The power-law is found to arise from quasi-discontinuous 
changes in the slope of the time series of the heat transport associated with the nucleation of dislocation 
pairs and roll pinch-off events. For larger frequencies, the power spectra decay exponentially as expected 
for time continuous deterministic dynamics. 
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Significant insight into the onset of chaotic dynamics in fluid systems, and continuum systems 
in general, has been gained from cryogenic Rayleigh-Benard convection experiments Q ||]. Two 
of the most dramatic discoveries were the observation of time dependence almost immediately 
above the onset of convective flow, and the power-law fall-off in frequency for the power spectral 
density derived from time series of a global measurement of the temperature difference across 
the fluid at fixed heat flow (3|, ||, S|. However, these and other important observations remain 
poorly understood. The power-law behavior was unexpected, since bounded deterministic models 
typically show an exponential falloff at high frequency |J. Phenomenological stochastic models 
were proposed to explain the spectra but no understanding of the origin of the ad hoc stochastic 
driving has followed. 

In this paper, we use numerical simulations of the three-dimensional Boussinesq equations for 
the fluid flow and heat transport in the cylindrical geometries of the experiments with realistic 
boundary conditions to investigate the power spectrum in more detail. The numerical simulations 
allow us to determine the spatial structure of the flow field in the aperiodic dynamics, and the 
absence of experimental or measurement noise provides us with more complete results for the power 
spectra. Our completely deterministic simulations yield results consistent with the experimental 
observations, including a power-law falloff of the power spectrum over the range accessible to the 
experiment. Using knowledge of the flow field, we are able to associate this power-law behavior with 
specific events in the dynamics, namely, the creation and annihilation of defects in the convection 
roll structure, which occur on a time scale rapid compared with the slow pattern evolution. At 
higher frequencies, the power spectra decay exponentially, consistent with the behavior expected 
for smooth deterministic time evolution. The low amplitude region of the spectra was inaccessible 
experimentally due to the noise floor. 

Our simulations in a cylindrical geometry are performed using an efficient spectral element 
algorithm (described in detail elsewhere ||). The velocity u, temperature T, and pressure p, 
evolve according to the Boussinesq equations, 

a' 1 (d t + u.v) u = -Vp + RTz + V 2 u, 
(d t + tf.v) T = V 2 T, 
V'U = 0, 

where dt indicates time differentiation, z is a unit vector in the vertical direction, a is the Prandtl 
number, and R is the Rayleigh number. The equations are nondimensionalized in the standard 
manner using the layer depth h, the vertical diffusion time for heat t v , and the constant temperature 
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difference across the layer AT, as the length, time, and temperature scales, respectively. All 
variables in the following discussion are nondimensional using this scaling. The lower and upper 
surfaces (z = 0, 1) are no-slip and are held at constant temperature. The sidewalls are no-slip 



and perfectly conducting [15], and the initial conditions are small random thermal perturbations 
of magnitude 0.2 imposed upon an otherwise quiescent layer, u = 0, T = 0. 

In nearly all cryogenic experiments, for reasons of increased experimental resolution, the heat 
flux across the convection layer, Q, and the temperature of either the upper or lower surface 
are held constant while measurements of AT(t) are made. These measurements are reported as 
R(t)/R c or AT(i)/AT c , where AT C is the temperature difference across the layer and R c is the 
Rayleigh number at the convective threshold. Theoretical calculations, on the other hand, most 
often consider both the upper and lower surfaces to be held at constant temperature and observe 
the time dependence in Q(t), which can be reported as a time series of the Nusselt number N(t) 
(the normalized heat current through the fluid layer). It has been shown experimentally that 
fixing Q or fixing AT does not appear to change the flow dynamics, and the conclusions from 
measurements of R(t) at fixed N or N(t) at fixed R will be similar ||. 

In order to make contact with experiment @, |9], |l(| we focus our discussion on simulations with 
aspect ratio T = 4.72 (T = r/h, r is the radius), a = 0.78 (experimental fluid was nonsuperfluid 
He 4 ), and constant AT. A key result of the experiments was the observation that the power 
spectrum, P{v), of measured R(t) values exhibited the power-law behavior, P(v) ~ v~ n ,n = 
—4.0 ± 0.2 over the frequency range 0.5 < v < 9 [JToJ] (results were reported for e = 3.62, where 
e = [R — R c )/R c is the reduced Rayleigh number). 

Six representative time series N(t) from our simulations are shown in Fig. [l|. In terms of the 
horizontal diffusion time for heat 7% (t^ = T 2 t v ), the simulation times are tf ~ lOOr^ (tf ~ 50^ 
for case (vi)), which is comparable with the longest experiments, tf ~ 65r/j (with one long run for 
tj 13577,) ||. This is considerably longer than IY^, which has been suggested as the earliest 
time scale for the flowfield to reach equilibrium [11|. However, as discussed below, we find that 
the dynamics can occur on even longer time scales. The simulated time-averaged values of N — 1 
are within 5.5% of the experimental values given by N - 1 = 1.034/3 + 0.981/3 3 - 0.866/3 5 , [3 = 
1 — R c /R Q. Spatial and temporal resolution studies have been performed to ensure the accuracy 
of the calculated values of N(t) for the chosen simulation parameters. 

We now consider the periodic time series in more detail, case (ii) in Fig. |l]. To determine the 
influence of the pattern dynamics on the power spectrum, we used a sliding window in time to 
calculate successive time-localized power spectra (a spectrogram) as shown in Fig. ||b. Most of 
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the power in the spectra can be attributed to nucleation of dislocation pairs (and, to a lesser 
extent dislocation annihilation). A plot of P{y) at a particular time (a vertical slice of Fig. |||b) 
yields the windowed power spectrum centered about that point in time in Fig. ||a. Figure || shows 
three such power spectra from the spectrogram evaluated with windows centered on (a), (d), and 
(e) corresponding to dislocation nucleation, annihilation, and glide, respectively. The local power 
spectrum centered on the nucleation of a dislocation pair generates a power-law region of significant 
magnitude; the local power spectrum centered on the dislocation annihilation generates a power- 
law region that is a factor of 10 smaller in magnitude; whereas the local power spectrum calculated 
during dislocation glide falls off more rapidly with v and does not make a significant contribution. 
The origin of the power-law is the time signature of the nucleation of a dislocation pair manifested 
in the quasi-discontinuous slope of N during the rapid excursion to a state of decreased heat 
transfer (an individual event is shown in Fig. |^). Note that a triangular feature with discontinuous 
changes in the slope yields a asymptotic behavior in the power spectrum. 

The power spectra of the remaining time series in Fig. [l] yield a power-law region P(y) ~ 
i/~ 4 because the dynamics (periodic, quasiperiodic, and chaotic) are dominated by the nucleation 
of dislocation pairs and roll pinch-off events. The differences between the chaotic and periodic 
dynamics are apparent in the low-frequency region and are more apparent on a log-linear plot of 
P(y). For example, Fig. || shows the similarity between N(t) for the periodic simulation, case (ii), 
and the chaotic simulation, case (vi). The dynamics in the chaotic state are much more complicated; 
however, they are dominated by the roll pinch-off events that maintain the characteristic quasi- 
discontinuous slope of N, yielding a power-law region in the power spectrum. A comparison of the 
power spectra for the periodic and chaotic time series is shown in Fig. [|. 

The average of the windowed power spectra of Fig. || eventually exhibit an exponential decay, 
which continues until reaching the spectrum floor; this is shown for case (ii) in Fig. |j[ Exponen- 
tial decay in the power spectra at high frequency is expected for bounded smooth deterministic 
dynamics . The exponential decay in the power spectra was not detected in experiment due to 
the presence of instrumental noise which masked the small scale region. 

In the cryogenic experiments, flow visualization was not possible leaving the precise details of 
the underlying pattern uncertain. With this in mind, we briefly discuss the dynamics represented 
in Fig. [|. Case (i) illustrates a time-independent Pan-Am pattern similar to panel (a) of Fig. ^. 
Case (ii) is periodic with period t = 8At\ (note the initial transient lasting 27r^); the dynamics of 
one period are illustrated in Figs. ^ and ||. Figure || displays the pattern at six different instances 
in time corresponding to the events labeled in Fig. ||a. Initially there is a Pan-Am pattern with 



4 



two opposing wall foci causing roll compression (a), eventually nucleating a dislocation pair in the 
center of the domain (b). The dislocations quickly climb to the wall (c), at which point they both 
begin to glide slowly toward the same wall focus. However, the lower dislocation is annihilated 
at the sidewall (d), and the remaining dislocation continues to glide slowly into the wall focus, 
where it is annihilated (e). A Pan- Am pattern again forms (f), and finally the process repeats. 
This is in general agreement with flow visualizations from related room-temperature argon-based 
experiments [12l p^| . Case (iii) may be periodic on a long time scale of t ~ 40^; the duration 
of the simulation is inadequate to be conclusive. Case (iv) illustrates a chaotic burst of duration 
t ~ 541^ bounded by periodic dynamics with a period of t ~ 17. Again the simulation duration 
is inadequate to determine whether this is a transient state or whether the chaotic bursting will 
repeat. Case (v) shows an initial chaotic transient that makes a transition at t ~ 18r^ to a very 
complicated quasiperiodic state where the central roll pair is pinned by the dynamic motion of 
two opposing disclinations. The dominant mode in the quasiperiodic state has a time scale t ~ 8. 
Case (vi) illustrates chaotic dynamics. 

We have also performed simulations for the T = 4.72 cylindrical domain with insulating lateral 
boundary conditions, in addition to simulations in a T = 7.66 domain (cr = 0.69 for argon) for both 
conducting and insulating lateral boundaries. Considering these additional results, we maintain 
our conclusions concerning the origin of the power-law. 

This work represents a joint computational and theoretical effort to further our quantitative 
understanding of complex dynamics in spatially extended nonequlibrium systems. An important 
link missing in nearly all theoretical work to date has been a quantitative comparison with experi- 
ment. Our results demonstrate that this quantitative comparison with experiment is now possible. 
We plan to use this approach to investigate spatiotemporal chaos in larger aspect ratio systems. 
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North Carolina Supercomputing Center. 
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FIG. 1: Plots of the dimensionless heat transport N(t) for cases (i-vi) for reduced Raylcigh number e = 
0.557,0.614,0.8, 1.0, 1.5, and 3.0, respectively. For cases (i-v), At = 0.01, and for case (vi), At = 0.005 {At 
is the time step). 
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FIG. 2: Time series N(t) (top) and corresponding spectrogram (bottom) for one period of case (ii). The 
labels a-f represent particular moments in the evolution of the pattern and are discussed in the text. The 
spectrogram displays 9 orders of magnitude of the power, P(v), with the smallest and largest contours 
labeled; the remaining contours each differ by a factor of 10. The spectrogram was calculated using a 
sliding Hann window of width At = 20.48 and linearly detrended overlapping segments (segments overlap 
by t = 20.0). 




FIG. 3: Windowed power spectra. The power spectra labeled (a), (d), and (e) are vertical slices of the 
spectrogram taken at representative times for case (ii) at dislocation nucleation, dislocation annihilation, 
and dislocation glide, respectively. The curves labeled (ii) and (vi) are the average of the windowed power 
spectra using the entire spectrogram for cases (ii) and (vi). The dashed lines represent P{v) ~ v~ A . 
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FIG. 4: A closeup of the time series N(t), illustrating the signature of a nucleation of a dislocation pair for 
case (ii) (solid line), and the signature of a roll pinch-off for case (vi) (dashed line). Time t* is measured 
as t* = (t — tj)/At; ti denotes when the event begins, and At is the duration of the event. For the events 
shown, ti = 587, 270.5 and At = 18, 1.5 for cases (ii) and (vi), respectively. 





FIG. 5: The power spectrum, (P(y)}, for case (ii) on a log-linear scale to illustrate region of exponential 
decay. The slope of the dashed line is -3.8; the crossover to exponential decay occurs at v w 1.5. 
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FIG. 6: Flow visualization showing contours of the thermal perturbation at the mid-depth, (6 evenly spaced 
contours; —0.2 < 6T < 0.2, negative values are dashed lines, and positive values are solid lines) for case (ii). 
Panels a-f are for t = 600, 605, 630, 650, 735, 785. The dislocations glide to the right; during the next period, 
the dislocations glide to the left, as can already be discerned in (f) by the bias in the roll compression. This 
left and right alternation continues for the entire simulation. 
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